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SUMMARY 

Friction laws, which are a key to the understanding of the diversity of earthquakes, are 
eonsidered theoretieally. Using dimensional analysis, the logarithmie dependenee of the 
frietion eoeffieient on the slip veloeity and the state variable is derived without any knowl¬ 
edge of the underlying physieal proeesses on the frietional surfaee. This is based on a 
simple assumption that the frietion eoeffieient is expressed as the differenee from a ref- 
erenee state. Therefore, the funetional form of the rate and state dependent frietion law 
itself does not neeessarily mean that thermal aetivation proeesses dominate frietion. It is 
also shown that, if there are two (or more) state variables having the same dimension, we 
need not assume the logarithmie dependenee on the state variables. 
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1 INTRODUCTION 

Numerous attempts have been made to clarify the frictional properties of rocks to understand the 
diverse nature of the rupture dynamics of earthquakes (iKanamori 19941) . They are investigated mainly 
using laboratory experiments, and some of them are summarized as an empirical law, which is now 
known as the rate and state dependent friction (RSF) law dPieterich 19781 IRuina 198^ . A general 
form of the RSF law describes the friction coefficient // as a function of the slip velocity V and the 
state variables {0} = 0i,02, ■ 




( 1 ) 
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The state variables represent any relevant microscopic features of the slipping interface. Convention¬ 
ally, the state variables are assumed to have the dimension of time. We thus assume that, throughout 
this paper, the state variables have the dimension of time unless otherwise indicated. 

Importantly, the RSF law involves the difference in friction coefficient between two states, and 
does not state anything about the absolute value of the friction coefficient. This is apparent in view of 
a concrete form of the RSF law (with a single state variable). 

iu(l/,6')= alog:^+61og^, (2) 

V if ^if 


where 14 and 0* are for a reference state, and a and b are nondimensional positive constants. 

Because the state variables {6} describe the microscopic properties on the slipping surface, they 
cannot be regarded as controllable parameters but time-dependent quantities. Thus, the RSF law in the 
form of Eq. © must be supplemented with time evolution laws for state variables. The time evolution 
of state variables depend on the microscopic physical processes on the frictional surface, and therefore 
the evolution laws involve the microscopic details of the system. In this paper, the time evolution laws 
are not discussed. 

In most cases, a single state variable version, Eq. ([21), has been successful in fitting experimental 
data on rock friction if the time evolution law for 9 is chosen suitably (IMarone 19981) . If an empirical 
law reproduces experiments well, it often gives a clue to the underlying (microscopic) physics. Ac¬ 
tually, Chester (1 19941 1 and Heslot et al. (1 19941 1 suggest that the RSE law is a consequence of thermal 
activation processes at asperities, and some experimental work followed to validate the proposition 
(INakatani 20011 ). In this respect, the logarithmic term with respect to the slip velocity in Eq. (O ap¬ 
pears natural. 

Here we show an alternative view on the RSE law; Eq. Q is derived using dimensional analysis. 
This means that the logarithmic term itself does not imply any microscopic physical processes. (Note 
that this does not mean that thermal activation processes is irrelevant to the RSE law.) In a similar 
manner, we also show the generally permissible form of the RSE law in the existence of two or more 
state variables. 

It should be stressed that the present attempt is not only of purely theoretical interest, but may con¬ 
tribute to the better understanding of RSE as well as to the formulation of new friction laws. The RSE 
law with multiple state variables has been used to explain some experimental data under hydrother¬ 
mal conditions (IGu et al. 19841 Blanpied & Tullis 1986} iNoda & Shimamoto 20121 ). but there is no a 
priori reason for assuming an equation like (O with multiple state variables. A clear demonstration 
of this result is given. In addition, recent numerous laboratory data show that the conventional RSE 


law should be replaced by a new law at seismic slip velocities (mm/s to m/s) (Goldsby & Tullis 2002 
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iHan et al. 2007l[Han et al. 2011llRice 2006I) . We discuss a possible application of the RSF formulation 
to this regime. 


2 ONE STATE-VARIABLE CASE 

Throughout this paper, we wish to focus on empirical friction laws. Hence, the friction coefficient must 
be expressed in terms of experimentally accessible quantities only. This assumption may be natural if 
one cannot access sufficient microscopic information for physical modeling. 

If the friction coefficient depends on the sliding velocity V and the state variable 0 that has the 
dimension of time, it can be formally written as 

p = p{v,e). (3) 

Here p{V,9) means the friction coefficient that is experimentally measured at a given velocity V and 
a certain value of the state variable 9. Although the state variable 9 itself is not directly measurable, 
one can instead compute 9 from other observable quantities using an evolution law. 


2.1 Some consequences of the H theorem 


Note that the friction coefficient is a nondimensional quantity. According to the H theorem for dimen¬ 
sional analysis (iBarenblatt 20031) . any nondimensional quantity must be a function of nondimensional 
parameters only. In view of Eq. ([3]), however, we cannot produce any nondimensional quantities us¬ 
ing velocity and time (i.e., the state variable). Therefore, the right hand side of Eq. (I3]l cannot be 
nondimensional unless other relevant parameters come into play. 

Eor example, one can introduce a length constant L that is defined as fhe characferisfic relaxation 
lengfh for friction. 


p = p{V,9,L). 


(4) 


Then V9/L is the only nondimensional quantity that can be constructed from V, 9, and L. Erom the 
H-theorem, one can write 


h = 



(5) 


where /(•) is a scaling function. Although Eq. ([51) is consistent with dimensional analysis, it cannot 
represent the slip velocity dependence of the steady-state friction coefficient. To see this, note that the 
state variable at steady state must be proportional to L/V, as there are no other constants that have the 
dimension of time. Then, obviously, V9/L is independent of the sliding velocity. 

To express both the steady-state slip velocity dependence and transient behavior, the two argu¬ 
ments in Eq. Q, V and 9, must be nondimensionalized by other dimensional parameters. Such pa- 
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rameters are usually material constants or other intrinsic parameters of a system. For example, the slip 
velocity V and the state variable 6 may be normalized using the sound velocity and the Debye 
frequency ujd, respectively. Then we can write instead of Eq. ([3]) 

F = / ■ (6) 

If this nondimensionalisation is appropriate, experimental data for many materials having different 
sound velocities and Debye frequencies should collapse to a single master curve. To the best of our 
knowledge, however, there have been no attempts to demonstrate such a dependence to this date. 
To choose an appropriate nondimensionalisation, we must know the essential physical processes be¬ 
hind a phenomenon in focus and, in general, must consider some models based on microscopic level 
physics. Such an approach would lead to an alternative form of the friction law without invoking a ref¬ 
erence velocity, which is conventionally denoted by 1/*. Although there are some promising attempts 
dChester 19941 IHeslot et al. 19941 IBar-Sinai et al. 20141) . this is beyond the scope of the present work. 
More detailed discussion on the microscopic nature of friction law will be given elsewhere. 


2.2 Nondimensionalisation with reference state 

2.2.1 General form 

If one cannot know intrinsic parameters that have the dimensions of velocity and time, nondimension¬ 
alisation of V and 9 can be achieved by introducing a reference state. Namely, one may use in place 
of intrinsic constants the slip velocity and the state variable at a reference state, 14 and 0*. Then we 
can utilise two nondimensional parameters: E/14 and 9/6^. The choice of a reference state {V^,9fj 
may be arbitrary, but it is convenient to choose 0* as the steady-state value at sliding velocity of E*. 
Then /r(E*, 0*) is just the friction coefficient measured at the steady state of sliding velocity E*. With 
the aid of a reference state, one can describe the friction coefficient in the form that is compatible with 
dimensional analysis. 

From the physical point of view, the friction coefficient at (E, 0) must not depend on (E*, 0*). One 
thus cannot write p = p{V, 9, E*, 9fj, but must distinguish p{V, 9) and ^(E*, 0*). Then a friction law 
should be written as 

F{p{V,9),p{V,,9,)) = f (7) 

where the functional form of is arbitrary. This is a general form of friction law that does not 

contain any intrinsic dimensional parameters. The choice of Ff, •) defines a concrete form of friction 
law. 
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2.2.2 Derivation of logarithmic dependence 


One may describe the friction coefficient as the change from a reference state value; i.e., F{x, y) = 
X — y. Then Eq. dV]) leads to 


p{v,9) -p{v,,e,) = f 


V e_ 
0* 


( 8 ) 


It can be shown that Eq. ([8]) is equivalent to the RSE law, Eq. Q- Differentiating the both sides of 
Eq. ([8]l with respect to V, one obtains 


dpiV,9) 


1 


d 




t'V_ 9_ 


(9) 


dV KWK) VK’0: 

As V, 14, 9, and 0* can be chosen arbitrarily, one may set = V and 9 = 0*. Then Eq. (|9l) becomes 


dp{V,9f) 1 df{si,S2) 


dV 


V dsi 


df{si,S 2 ) 


dsi 


a 

V' 

so 


( 10 ) 

( 11 ) 


where the scaling function / is denoted by /(si, S 2 ) and sq = (1,1). Noting that a is a constant, one 
can integrate Eq. (ITOl) with respect to V from 14 to V and obtain 


V 

p.{V-,9f) -/i(14,('*) = olog—. 

V if: 

In the same manner, the following equation can be easily derived. 


p{V,9) - p{V,9= 51og—, 

df{si,S2) 


ds2 


( 12 ) 

(13) 

(14) 


so 


Adding Eq. (fT^ to Eq. ([T3]) leads to Eq. (Illl. Namely, Eq. ([8]l is equivalent to Eq. (O. 


2.2.2 Another example 

Note that Eq. ([8]l is not the only consequence of the dimensional analysis. Alternatively, one may 
choose the ratio of friction coefficients measured at two different states; i.e., F{x, y) = xjy. Then Eq. 
(ITI) leads to 



Erom this equation, one can easily derive 


p{V,9) = piV,,9,) (16) 

where the exponents a and /3 are constants. 
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2.3 Summary of one state-variable case 

A choice of F(-, •) in Eq. ([71) is arbitrary, but it determines the functional form of friction law in¬ 
evitably. If we choose to describe the friction coefficient as the difference from the reference state, the 
logarithmic form of friction law is inevitable as a direct consequence of dimensional analysis. This 
constitutes our central result. Other choices are also possible, but they yield different forms for the 
friction law. 


3 THE CASE OF TWO STATE VARIABLES 

The RSF law is sometimes extended to include a second state variable. Hereafter the second state 
variable is denoted by cf). We assume that the friction coefficient is uniquely determined by these three 
parameters. Namely, we can write 

fi = ( 17 ) 

The possible form of the friction law depends entirely on the dimension of the second state variable, 
(j). (Here we assume that 9 has the dimension of time.) We discuss two cases separately. 


3.1 If-(j) has the dimension of time 


First we discuss the case in which cf) has the dimension of time. As the dimension of 9 is also time, the 
ratio 0/(/) is nondimensional. Then the H theorem concludes /r = f{9/(l>), where / is an arbitrary func¬ 
tion. In this case, the evolution laws for 9 and 0 must be chosen appropriately so as to fit experimental 
data. However, we are unaware of such a formulation of friction law. 

Alternatively, one may again introduce a reference state, which is denoted by (K, 0*, (/>*). In this 
case, there are actually four independent nondimensional numbers: V/K> 9/9^,, ^/(/>*, 9/(\). Then 
following the same logic as in the previous section, one is led to 

V 9 6 


^{V, 9,4>) - /x(K, 6**, (/>*) = / 


(18) 


The scaling function / on the right hand side will be denoted by /(s) or /(si, • • •, S 4 ). Then, differ¬ 
entiating the both sides of Eq. (IT^ with respect to 9, and letting V = V^, 9 = 9^, and cj) = 4>*^ on® 
obtains 


dfJ-jV, 9, (j)) 
89 


1 ^ 
9 ds2 


so 


1 ^ 

cj) dS4^ 


(19) 


so 


where sq = {l,l,l,9/4>). Note that the coefficients df /ds2\so and df /ds 4 \so depend on 9/(j). There¬ 
fore, one cannot conclude 9, (f>) — 9^,(j)) oc log(0/0*) from this equation. 

One can repeat the same discussion with respect to (j) and has 
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dp{V,OA) 

d(j) 


1 ^ 
4> ds3 


so 


A K 

(jF dsi 


so 


Again, this equation does not lead to p,{V, 9, f) — p{V, 9, ff) oc log((^/(/)*). 
Contrastingly, the following equation can be obtained with respect to V. 
dp(V,9,(/>) 1 df 

dv V dsi so ’ 

and therefore 


fi{V, 9, f) - p{V^,9,4>) = a log y, 




K 

dsi 


so 


Taking Eqs. (fT^ . (l20l) . and (|2^ into account, we obtain 


p{V, 9, f) - p{V^,9^,4>^) = a 




^ A ^ 

9f ff (p 


( 20 ) 


( 21 ) 


( 22 ) 

(23) 


(24) 


where G is an arbitrary function. 

Two state variable RSF law is often assumed as the following form (IGu et al. 19841 Blanpied & Tullis 1986 
INoda & Shimamoto 2012I) . 


p{V,9,(j)) - p{V^,9^,4>^) = alog-^ + 6 ilog^ + 62 log (25) 

V* 

Note that this is just a special case of Eq. (|24]) . One may thus assume Eq. (|25] ) as long as it can 
describe experimental data well, because the ultimate justification of friction laws relies entirely on 
the comparison with experiments. Here we just point out that Eq. (|2^ is more general than Eq. (|25] ) 
from the viewpoint of dimensional analysis. 


3.2 If the dimension of f is independent of V and 9 


We assume that there exist no set of integers (n, m, 1) that V^9'^(jP is nondimensional. If such integers 
exist, the H theorem concludes p = f{V'^9'^(lF), and no restriction is assigned to its functional form 
from dimensional analysis only. 

If such integers do not exist, one may write the analogue of Eq. Q from dimensional consistency. 


p{V,9,(t)) - p{V^,9^,(j)^) 


f 


/V 9 f 


(26) 


Following the same discussion as the one state variable case, the above equation leads to Eq. (1251) . 
Therefore, Eq. (1251) is a general form of friction law if one cannot make a nondimensional number 
using V, 9, and f. 

An example of such f is the temperature, the stresses, or the mass density. Particularly, adopt¬ 
ing the temperature as the second state variable may be plausible in modeling phenomena in which 
frictional heating is relevant. 
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4 THREE OR MORE STATE VARIABLES 

The state variables are denoted by Then, if there exist no such integers (n, m*) that H ^1"* 
nondimensional, and if no intrinsic physical constants are known, the deviation of the friction coeffi¬ 
cient from a reference state, /j,(V* , {0 *}), must be written as 

f,(V, {0,}) - {0*}) = «log + E log (27) 

However, such choice of state variables are generally impossible if the number of state variables is 
larger than three, as the number of independent dimensions are typically four: mass, time, length, 
and temperature. If there are four state variables, all the variables including the velocity cannot be 
independent and therefore Eq. (|27]) does not necessarily hold. Then, as in section [TTI no restriction is 
imposed to the functional form of friction law from dimensional analysis alone. 


5 CONCLUSIONS 


In this paper, with the aid of dimensional analysis, it is shown that the logarithmic terms in the RSF 
law are inevitable irrespective of the actual physical processes if one describes the friction coefficient 
in terms of the difference from the reference-state value. With regard to the two state-variable formu¬ 
lation, the logarithmic terms must be interpreted as an additional assumption unless two state variables 
have independent dimensions. 

In general, the formulation of empirical laws with a reference state is rather natural and univer¬ 
sal when one cannot know the relevant physical processes. Therefore, the formulation presented here 
may be still useful in describing novel frictional properties, for example, significant weakening be¬ 
haviors at intermediate to high slip velocities, mm/s to m/s ( Goldsby & Tullis 2002} [Han et al. 20071 
IHan et al. 2011tldi Toro et al. 201111 . In these phenomena, the second state variable may represent a lo¬ 
cal temperature at microscopic junctions (iRice 20061) . the mass density of silica gel (Goldsby & Tullis 2002 1 , 
or the density of nanopowders (IHan et al. 20111) . In each case, the dimension of a second state variable 
is independent of time. Therefore, the discussion given in section [T2] applies and we may assume Eq. 

(1251) as the empirical friction law. 

We conclude that the logarithmic term is derived if (1) one does not know an intrinsic velocity con¬ 
stant and uses a reference state velocity, and if (2) the friction coefficient is expressed as the difference 
from a reference state. Then, the contraposition is the following: A non-logarithmic behaviour implies 
that (3) an intrinsic velocity constant is used and a reference state velocity is ruled out, or (4) the fric¬ 
tion coefficient is not expressed as the difference from a reference-state value. For instance, a velocity¬ 
weakening behaviour in some experiments at high slip velocities (Idi Toro et al. 201 llIKuwano et al. 20TTI) 
might be described as 
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pi{V) = po- (28) 

where po, A/i, and are constants (iRice 20061) . This equation contains the intrinsic velocity constant 
14), and therefore corresponds to case (3) above. An example for case (4) is Eqs. lfTSl) and (fT^ . where 
the friction coefficient is expressed as the ratio to a reference-state value. 
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